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Abstract 

The concepts of quantile position, trajectory, and velocity are denned. For a tunnel- 
ing quantum mechanical wave packet, it is proved that its quantile position always 
stays behind that of a free wave packet with the same initial parameters. In quan- 
tum mechanics the quantile trajectories are mathematically identical to Bohm's 
trajectories. A generalization to three dimensions is given. 
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The discussion of tunneling times has a long history in the literature, for a 
review see [1]. Recent developments are, e.g., the "tunneled flux" approach [2], 
the operational projector approach to tunneling times [3], and the calculation 
of tunneling times in the framework of the Dirac equation [4]. Among other 
findings it has been reported that the velocity of a particle can be larger within 
a repulsive barrier than outside the barrier and that instantaneous tunneling 
can occur. Great care has to be exercised, however, in denning a velocity in 
quantum mechanics, since the velocity definition of classical mechanics re- 
quires the concept of a trajectory for a point particle which breaks down in 
quantum mechanics. In this letter we introduce a definition of a velocity which 
is strictly based on probability concepts. 
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1 Definitions 



For any probability density g(x) the quantile xq associated with the proba- 
bility Q is defined in the mathematical literature, c.f., e.g., [5,6], by 

X Q 

Q = J q{x) dx . 



For a position- and time-dependent probability density g(x, t) we introduce 
the time-dependent quantile position xp(t) through 

g(x, t) dx . (1) 

X P (t) 



The change to the complementary integration interval is chosen merely for 
convenience. Essential is the transfer of the quantile concept from a probability 
density in mathematical statistics to a space- and time-dependent distribution 
describing a physics problem. It yields for every time t a well-defined point 
xp(t). Considering x = xp(t) as a quantile trajectory we define the quantile 
velocity vp(t) = dxp/dt, which obviously depends on the chosen value of P. 

In the case of a conserved probability the corresponding probability density 
fulfills the continuity equation of the form 

g<«, ( ) + !<*.«) = 0. (2) 



j(x,t) being the corresponding probability current density vanishing at infin- 
ity. Equations (1) and (2) permit a direct calculation of the quantile velocity 
vp(t). Differentiating (1) with respect to time we immediately obtain 

dx P j(x P ,t) 

v p — — — —, r • (<j) 

dt g(x P ,t) v ' 



Equation (3) is an ordinary differential equation for the quantile position xp(t), 
an implicit solution of which is given by (1). 

In an experiment with quantum-mechanical wave packets the quantile velocity 
can be determined on a statistical basis by time-of-flight measurements: One 
prepares by the same procedure N single-particle wave packets and sets a 
clock to zero at the moment at which the spatial expectation value of a wave 
packet leaves the source. With a detector placed at position x x one registers 
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the arrival times t\ m of particles for m = 1,2, ... ,N and orders them such 
that tn < t±2 < ■ ■ •• One picks the time t\ n which is the largest of the smallest 
n times and chooses n/N = P. The time t\ n is the arrival time of the quantile 
xp at the position x\, i.e., xp(ti n ) — x±. By repeating the experiment with a 
detector at x 2 one obtains t 2n , etc. The points xp(t in ) are discrete points on 
the quantile trajectory x — xp(t). If x\ and x 2 mark the beginning and the 
end of a potential barrier then t 2n — t\ n is the quantile traversal time of the 
barrier. 

In the case of classical systems, e.g., classical electrodynamics, the quantile 
velocity can quite naturally be interpreted as velocity of an electromagnetic 
signal. We identify the density g(x,t) = w(x,t)/W with the ratio of the elec- 
tromagnetic energy density w and the total energy W of the pulse. We say the 
signal has left the transmitter, if the fraction P of the total energy has left the 
transmitter. It has reached the detector if the fraction P has been absorbed by 
the detector. Obviously a minimum amount of energy (the threshold energy 
Wtb) is needed in the detector to register a signal. One will therefore choose 
P such that PW > W th . 



2 Free wave packet 

As an example we consider the time development of a time-dependent Gaus- 
sian probability distribution 

This probability density describes just as well (i) the marginal distribution of a 
bivariate Gaussian phase-space distribution of a classical assembly (see, e.g., 
[7]) of force-free particles of mass m with initial position expectation value 
x, initial spatial variation a xQ , and momentum expectation value p = mv, 
momentum variation a p = h/(2<T x0 ), and velocity variation a v = a p /m, and 
(ii) the spatial probability distribution of a quantum-mechanical force-free 
Gaussian wave packet with the same initial parameters as the above classical 
distribution. Equation (1) yields the quantile trajectories 

x P (t) = x + vt + (<T x (t)/a x0 )(x - x) , 

and thus the quantile velocities 

vp(t) =v + (a 2 v t/(o x (t)o x0 )) (x -x) . 
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Here xq = xp(0) is the initial quantile position associated with P according to 
the condition (1/2) erfc[(xo — x)/(y/2a x o)] = P. If we take (4) to be a classical 
phase-space distribution, we note that the quantile trajectories x = xp{t) 
are not trajectories of free particles possessing constant velocities v = p/m. 
Only for P = |, i.e., Xq = x is the quantile trajectory identical to a particle 
trajectory. Quantile trajectories of a free wave packet for different values of P 
are shown as dotted lines in Fig. 1. 



3 Quantile motion for a non-conserved probability 



Our definition (1) allows a consistent description of motion also in the case 
of a non-conserved probability. We assume that instead of a normalization to 
one at all times, we have a loss of probability, 

oo 

J g(x,t)dx = F(t) < 1 , F(0) = 1. 

— oo 



As an illustration we consider again the time development of a wave packet 
prepared at t = to be a Gaussian distribution with the same parameters 
as in (4). In addition we introduce a temporal exponential decrease of the 
probability to find the particle anywhere in space, as for example caused by 
the presence of a constant purely imaginary potential. Thus the probability 
density may be written in the form 

I \ 1 f (x — X — Vt) 2 \ r i t /r\ 



where A is the temporal rate of probability loss and a x (t) is given in equation 
(4). Instead of the continuity equation (2) we have now 

^(x,t) + ^(x,t) = -l(x,t), (6) 



jx(x, t) being the current density of the non-conserved probability and l(x, t) = 
\q\(x, t) > describing the density of the temporal loss rate. For the quantile 
position xp(t) we obtain the integro-differential equation 



dx P (t) jx(x P ,t) 



dt Q\(x P ,t) Q\(x P ,t) 

Xp(t) 



l(x',t)dx'. (7) 
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The solid lines in Fig. 1 are quantile trajectories for different values of P. The 
curves end at the quantile position xp max where the integral over the loss rate 
on the right-hand side of (7) equals the probability current density j\(x,t). 
At later times the probability loss term in (7) dominates over the probability 
current density j\(x,t). For positions at values x > xp mSuX the condition (1) 
can no longer be satisfied. Note that the formulation of quantile motion for a 
non-conserved probability has general validity, since Eqs. (6) and (7) do not 
depend on the details of the probability loss. 



4 Tunneling 



In the following we show that for the same value of P and for two identi- 
cally prepared, tunneling and free wave packets, the quantile position of the 
tunneling wave packet in the transmission region of the potential remains at 
any given time t behind the quantile position of the free wave packet, i.e., 
the arrival time of the tunneled quantile behind the barrier is always later 
than that of the free wave packet. We consider the tunneling of a wave packet 
with spectral function ip{k) in momentum space through a repulsive square 
potential barrier, V(x) = V > for \x\ < a and V(x) =0 for \x\ > a. The 
probability densities of the tunneling and of the free wave packet with the 
same spectral function ip{k) are denoted by Qt(x, t) and g F (x,t), respectively. 
The functions 

oo oo 

P T (x,t) = J g T (x',t)dx' , P F (x,t) = J g F (x',t)dx' 

X X 



give the probabilities at time t for the tunneling and the free particle to be 
found to the right of the position x in the transmission region x > a of 
the potential. Since Pt(x, t) and Pp(x,t) decrease monotonically with x, our 
statement is always true if the condition 

AP(x,t) = P F (x,t) - P T (x,t) > (8) 



for x > a is fulfilled. We consider a wave packet described by a spectral 
function ip(k) with a positive wave number spectrum, i.e., ipik) =0 for k < 0. 
The difference AP(x, t) in the region x > a is 



AP(x,t) = ±J /[I - T(k')T*(k")] ^ dk ' dk " ' 

oo 
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where ip(k; x, t) = ip(k) exp[ik(x — x) — ihkH / (2m)} and 



T{k) = 



4k<y 



V(k) — (k + 7 ) 2 e 2ia(/c ' 7) — (k — 7 ) 2 e 2ia(fc+7) , 



is the transmission amplitude in momentum space for the repulsive square 
potential barrier (7 = [k 2 — 2mV /U 2 ) 1 ! 2 ). We consider AP(x,t) as the limiting 
function of a series of parameterized functions 

Ap( T t)-±7 /^MlW 
Al , J "2fi J V(k')V*(k") 

00 v ' v ' 

x $(k';x,t)r(h";x,t) dk , dk „ 
k' - k" + is 



for A — > 1, where P A (jfe) = (A; + 7 ) 2 e 2iaA ( fc ^ - (A; - 7 ) 2 e 2iaA ( fe+ ^. Thus, we have 
AP\(x,t) = for A = 0, so that the difference of probabilities is 



AP(x,t) = J—AP x (x,t) dA. 



(9) 



Explicitly calculating the derivative in (9) yields 



dA 



+4 (— ) 



{AmV\ 





00 



3 2iaAfc s i n (2aA 7 ) 



2?(fc) 



e 2iafe sin (2 a7 ) 

W) 



ip(k] x, t) dk 



dA 



ip(k; x', t) dk 



dx' 



(10) 



where Z x (k) = (A; + 7 )e 2iaA ( fc -^ - (A;- 7 )e 2iaA ( fc+ T). The positivity of (10) proves 
the condition (8) and therefore the initial statement on the retardation of the 
quantile trajectory of the tunneled wave packet in the transmission region 
of the potential. This procedure can be extended to the case of a general 
non- negative potential V(x) of finite range, introducing a segmentation of the 
potential into thin square potential barriers and then similarly parameterizing 
the difference AP(x,t) to obtain a positive-definite expression [8]. 

In Fig. 2a the time development of a wave packet incident on a barrier is 
shown. The packet is partly reflected and partly transmitted. In Fig. 2b the 
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quantile trajectories x = xp(t) for the same problem are presented. For small 
values of P the trajectories penetrate the barrier, for larger values they reverse 
their direction. We observe that in all cases the quantile velocity within the 
barrier region is smaller in absolute value than far from the barrier. 



5 Relation to Bohmian mechanics 

Bohm's interpretation of quantum mechanics [9] introduces particle trajecto- 
ries x — x(t) satisfying the equation of motion dx(t)/dt = j(x(t),t)/g(x(t),t). 
For the solution of this differential equation the initial condition x(0) = Xq 
is needed as a hidden parameter. If we set Bohm's initial position xq equal 
to the initial quantile position xp(0) the quantile trajectories are mathemat- 
ically identical to Bohm's particle trajectories. Conceptually, however, they 
are based on the probability interpretation of the conventional quantum me- 
chanics. 

Leavens and Aers [10] and Spiller et al. [11] realized that Bohm's interpretation 
of quantum mechanics — reintroducing the point-particle concept — offers the 
opportunity of unambiguously defining the velocity as in classical mechanics 
(see also [12,13]). In contrast to their work our approach is based on the 
conventional probabilistic concepts and does not depend on the interpretation 
given by Bohm. In Ref. [10] tunneling times are given which are based on the 
concept of Bohm trajectories. McKinnon and Leavens [14] pointed out the 
significance with respect to tunnel times of a particular Bohm trajectory for 
P = |T| 2 , with |T| 2 being the transmission probability. Wu and Sprung [15] 
have discussed quantum probability patterns for stationary problems. 



6 Generalization to three dimensions 

The equivalence of equations (1) and (3) has been noticed before, see, e.g., 
[12,16]. In [12] it was considered as an artifact of low dimensionality. In both 
papers, this equivalence was not related to the probability interpretation made 
possible by equation (1). However, the generalization to three dimensions is 
straightforward, if one considers that the interval between two quantile posi- 
tions xp^t) and Xp 2 (t) always contains the same amount of probability Q = 
Pi —P 2 and that all trajectories with initial positions x Pl (0) < x P (0) < xp 2 (0) 
stay in the interval x Pl (t) < x P (t) < x P2 (t). This suggests that in three di- 
mensions trajectories defined by the equation of motion 




(11) 
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and having initial positions xp(0,x) = x within a volume V containing an 
amount P of probability stay within a volume V t containing the same amount 
of probability. To prove this statement we introduce the time- dependent base 
of tangent vectors rf(t,x) = dxp(t : x)/dxi, where % = 1,2,3. Their Jacobian 
is given by the expression D(t : x) = r7 1 (t,x) • [rj 2 (t,x.) x rj 3 (t, x)], whereas its 
time derivative satisfies 

3D 

— (t,x) = £>(f,x)V Xp • v(x P (f,x),*) . 

In particular at the initial time t = we have D(0) = 1. The total time 
derivative of the density is then given by 

dg(x P ,t) 

= -g(x P , t) V Xp • v(x P , t) , 



so that we find the product p(t, x) = g(x. P (t, x), t)D(t, x) to be time indepen- 
dent, 

|«,x) = 0. 

This implies that an integral over a finite volume Vo, 

| pit, x') d 3 x' = P , 
Vo 

is time independent. An amount of probability 

P = J p(0,x')d 3 x' = y^(x,0)d 3 x 

Vo V) 

contained initially in a volume V is at time t given by 

P = J Q(x P (t,x.'),t)D(t,x. f ) d 3 x. 
V) 

The time-dependent substitution x = xp(t,x') yields 



P = J e ( x , t) d 3 x 
v 



with the volume given by all the points xp(t, x) at time t lying on trajecto- 
ries with initial positions xp(0, x) = x in the volume Vq. For a simple example 
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quantile trajectories of a three-dimensional wave packet are shown in Fig. 3. 
Requiring the quantity P to be time independent, is a necessary condition for 
the validity of (11). Sufficient conditions will be discussed in a forthcoming 
publication. We note that a corresponding statement holds also for a system 
with a loss term in the continuity equation, c.f. (6). 



7 Concluding remarks 

We wish again to emphasize that the concept of quantile velocities is not lim- 
ited to problems of classical statistical assemblies and problems of quantum 
mechanics but can easily be extended to describe the propagation of pulses of 
electromagnetic energy in vacuum, in dispersive and absorptive media, and in 
wave guides [8]. Trajectories derived from (11) for electromagnetic phenomena 
(with j and g being the energy current density and the energy density, respec- 
tively) were already discussed by Holland [17] but no quantile interpretation 
was given. 
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Fig. 1. Quantile trajectories for the wave packet (4) (dotted lines) and (5) (solid 
lines). Trajectories are plotted for values of P between 0.1 and 0.9 in steps of 
AP = 0.2. In all figures positions are given in units h/y/eVm and times in units 
h/eV with m being the mass of the particle. At time t = the wave packet is in 
both Gaussian wave packet with mean position x = — 10h/y/eVm, mean 

momentum p = 2\J eVm, and momentum width a p = 0.2\/eVm. The wave packet 
(5) is additionally characterized by a probability loss of temporal rate A = O.lh/eV. 
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Fig. 2. (a) Time development of a Gaussian wave packet with the same initial pa- 
rameters as in Fig. 1 incident onto a repulsive potential barrier (indicated as broken 
line) of height Vq = 10 eV extending from x = —0.3h/y/eVm to x = 0.3h/y/eVm, 
where m is the particle mass. The individual graphs show the probability density 
q(x, t) as function of x for fixed values of t. The hatched area comprises the prob- 
ability P = 0.25. The thick line is the quantile trajectory for that value of P. (b) 
Quantile trajectories x = xp(t) (solid lines) for values between P = 0.1 and P = 0.7 
in steps of 0.05 for the wave packet shown in (a) but presented in a smaller part 
of the x,t plane. The thicker line corresponds to the trajectory shown in (a). The 
barrier region is bounded by broken lines. For comparison the quantile trajecto- 
ries x = Xp Tec \t) of a free wave packet are also shown as dotted lines. Note that 
x P (t) < x { p cc \t) for every pair of values P and t. 
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Fig. 3. The expectation value of a free three-dimensional spherically symmetric 
Gaussian wave packet, which initially (at time t = 0) is at the origin, moves in 
the positive x direction. As initial volume Vj comprising the probability P a sphere 
around the origin is chosen. Quantile trajectories xp(i) of points which at t = 
lie on the surface of Vo at later times lie on the surface of a volume Vt which 
also comprises the same probability P. In this simple example all volumes Vt are 
spheres. The plot shows the cuts z = through two spheres Vq, Vt which are circles 
and trajectories in the x, y plane. 
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